Why Are Left-Handed G-Quadruplexes Scarce?

G-quadruplexes (G4s) are nucleic acid structures crucial for the regulation of gene expression and genome maintenance. While they hold promise as nanodevice components, achieving desired G4 folds requires understanding the interplay between stability and structural properties, like helicity. Although right-handed G4 structures dominate the experimental data, the molecular basis for this preference over left-handed helicity is unclear. To address this, we employ all-atom molecular dynamics simulations and quantum chemical methods. Our results reveal that right-handed G4s exhibit greater thermodynamic and kinetic stability as a result of favorable sugar-phosphate backbone conformations in guanine tracts. Moreover, while hydrogen-bonding patterns influence helicity-specific G4 loop conformations, they minimally affect stability differences. We also elucidate the strong correlation between helicity and the strand progression direction, essential for G4 structures. These findings deepen our understanding of G4s, providing molecular-level insights into their structural and energetic preferences, which could inform the design of novel nanodevices.

How acidic amino acid residues facilitate DNA target site selection.

Despite the negative charge of the DNA backbone, acidic residues (Asp/Glu) commonly participate in the base readout, with a strong preference for cytosine. In fact, in the solved DNA/protein structures, cytosine is recognized almost exclusively by Asp/Glu through a direct hydrogen bond, while at the same time, adenine, regardless of its amino group, shows no propensity for Asp/Glu. Here, we analyzed the contribution of Asp/Glu to sequence-specific DNA binding using classical and ab initio simulations of selected transcription factors and found that it is governed by a fine balance between the repulsion from backbone phosphates and attractive interactions with cytosine. Specifically, Asp/Glu lower the affinity for noncytosine sites and thus act as negative selectors preventing off-target binding. At cytosine-containing sites, the favorable contribution does not merely rely on the formation of a single H-bond but usually requires the presence of positive potential generated by multiple cytosines, consistently with the observed excess of cytosine in the target sites. Finally, we show that the preference of Asp/Glu for cytosine over adenine is a result of the repulsion from the adenine imidazole ring and a tendency of purine-purine dinucleotides to adopt the BII conformation.

Prediction of protein assemblies, the next frontier: The CASP14-CAPRI experiment.

We present the results for CAPRI Round 50, the fourth joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of twelve targets, including six dimers, three trimers, and three higher-order oligomers. Four of these were easy targets, for which good structural templates were available either for the full assembly, or for the main interfaces (of the higher-order oligomers). Eight were difficult targets for which only distantly related templates were found for the individual subunits. Twenty-five CAPRI groups including eight automatic servers submitted ~1250 models per target. Twenty groups including six servers participated in the CAPRI scoring challenge submitted ~190 models per target. The accuracy of the predicted models was evaluated using the classical CAPRI criteria. The prediction performance was measured by a weighted scoring scheme that takes into account the number of models of acceptable quality or higher submitted by each group as part of their five top-ranking models. Compared to the previous CASP-CAPRI challenge, top performing groups submitted such models for a larger fraction (70-75%) of the targets in this Round, but fewer of these models were of high accuracy. Scorer groups achieved stronger performance with more groups submitting correct models for 70-80% of the targets or achieving high accuracy predictions. Servers performed less well in general, except for the MDOCKPP and LZERD servers, who performed on par with human groups. In addition to these results, major advances in methodology are discussed, providing an informative overview of where the prediction of protein assemblies currently stands.

Modeling protein structures with the coarse-grained UNRES force field in the CASP14 experiment.

The UNited RESidue (UNRES) force field was tested in the 14th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP14), in which larger oligomeric and multimeric targets were present compared to previous editions. Three prediction modes were tested (i) ab initio (the UNRES group), (ii) contact-assisted (the UNRES-contact group), and (iii) template-assisted (the UNRES-template group). For most of the targets, the contact restraints were derived from the server models top-ranked by the DeepQA method, while the DNCON2 method was used for 11 targets. Our consensus-fragment procedure was used to run template-assisted predictions. Each group also processed the Nuclear Magnetic Resonance (NMR)- and Small Angle X-Ray Scattering (SAXS)-data assisted targets. The average Global Distance Test Total Score (GDT_TS) of the 'Model 1' predictions were 29.17, 39.32, and 56.37 for the UNRES, UNRES-contact, and UNRES-template predictions, respectively, increasing by 0.53, 2.24, and 3.76, respectively, compared to CASP13. It was also found that the GDT_TS of the UNRES models obtained in ab initio mode and in the contact-assisted mode decreases with the square root of chain length, while the exponent in this relationship is 0.20 for the UNRES-template group models and 0.11 for the best performing AlphaFold2 models, which suggests that incorporation of database information, which stems from protein evolution, brings in long-range correlations, thus enabling the correction of force-field inaccuracies.

Pseudopotentials for coarse-grained cross-link-assisted modeling of protein structures.

Pseudopotentials for the chemical cross-links comprising the glutamic- and aspartic-acid side chains bridged with adipic- (ADH) or pimelic-acid hydrazide (PDH), and the lysine side chains bridged with glutaric (BS2 G) or suberic acid (BS3 ) for coarse-grained cross-link-assisted simulations were determined by canonical molecular dynamics with the Amber14sb force field. The potentials depend on the distance between side-chain ends and on side-chain orientation, this preventing from making cross-link contacts across the globule in simulations. The potentials were implemented in the UNRES coarse-grained force field and their effect on the quality of models was assessed with 11 monomeric and 1 dimeric proteins, using synthetic or experimental cross-link data. Simulations with the new potentials resulted in improvement of the generated models compared to unrestrained simulations in more instances compared to those with the statistical potentials.

Mechanism of recognition of parallel G-quadruplexes by DEAH/RHAU helicase DHX36 explored by molecular dynamics simulations.

Because of high stability and slow unfolding rates of G-quadruplexes (G4), cells have evolved specialized helicases that disrupt these non-canonical DNA and RNA structures in an ATP-dependent manner. One example is DHX36, a DEAH-box helicase, which participates in gene expression and replication by recognizing and unwinding parallel G4s. Here, we studied the molecular basis for the high affinity and specificity of DHX36 for parallel-type G4s using all-atom molecular dynamics simulations. By computing binding free energies, we found that the two main G4-interacting subdomains of DHX36, DSM and OB, separately exhibit high G4 affinity but they act cooperatively to recognize two distinctive features of parallel G4s: the exposed planar face of a guanine tetrad and the unique backbone conformation of a continuous guanine tract, respectively. Our results also show that DSM-mediated interactions are the main contributor to the binding free energy and rely on making extensive van der Waals contacts between the GXXXG motifs and hydrophobic residues of DSM and a flat guanine plane. Accordingly, the sterically more accessible 5'-G-tetrad allows for more favorable van der Waals and hydrophobic interactions which leads to the preferential binding of DSM to the 5'-side. In contrast to DSM, OB binds to G4 mostly through polar interactions by flexibly adapting to the 5'-terminal guanine tract to form a number of strong hydrogen bonds with the backbone phosphate groups. We also identified a third DHX36/G4 interaction site formed by the flexible loop missing in the crystal structure.

The Product of Matrix Metalloproteinase Cleavage of Doxorubicin Conjugate for Anticancer Drug Delivery: Calorimetric, Spectroscopic, and Molecular Dynamics Studies on Peptide-Doxorubicin Binding to DNA.

Matrix metalloproteinases (MMPs) are extracellular matrix degradation factors, promoting cancer progression. Hence, they could provide an enzyme-assisted delivery of doxorubicin (DOX) in cancer treatment. In the current study, the intercalation process of DOX and tetrapeptide-DOX, the product of the MMPs' cleavage of carrier-linked DOX, into dsDNA was investigated using stationary and time-resolved fluorescence spectroscopy, UV-Vis spectrophotometry and isothermal titration calorimetry (ITC). The molecular dynamics (MD) simulations on the same tetrapeptide-DOX…DNA and DOX…DNA systems were also performed. The undertaken studies indicate that DOX and tetrapeptide-DOX can effectively bond with dsDNA through the intercalation mode; however, tetrapeptide-DOX forms less stable complexes than free DOX. Moreover, the obtained results demonstrate that the differences in DNA affinity of both forms of DOX can be attributed to different intercalation modes. Tetrapeptide-DOX shows a preference to intercalate into DNA through the major groove, whereas DOX does it through the minor one. In summary, we can conclude that the tetrapeptide-DOX intercalation to DNA is significant and that even the lack of non-specific proteases releasing DOX from the tetrapeptide conjugate, the presence of which is suggested by the literature for the efficient release of DOX, should not prevent the cytostatic action of the anthracycline.

Why do G-quadruplexes dimerize through the 5'-ends? Driving forces for G4 DNA dimerization examined in atomic detail.

G-quadruplexes (G4) are secondary structures formed by guanine-rich nucleic acid sequences and shown to exist in living cells where they participate in regulation of gene expression and chromosome maintenance. G-quadruplexes with solvent-exposed guanine tetrads show the tendency to associate together through cofacial stacking, which may be important for packaging of G4-forming sequences and allows for the design of higher-order G4 DNA structures. To understand the molecular driving forces for G4 association, here, we study the binding interaction between two parallel-stranded G-quadruplexes using all-atom molecular dynamics simulations. The predicted dimerization free energies show that direct binding through the 5'-G-tetrads is the most preferred of all possible end-to-end stacking orientations, consistently with all available experimental data. Decomposition of dimerization enthalpies in combination with simulations at varying ionic strength further indicate that the observed orientational preferences arise from a fine balance between the electrostatic repulsion of the sugar-phosphate backbones and favorable counterion binding at the dimeric interface. We also demonstrate how these molecular-scale findings can be used to devise means of controlling G4 dimerization equilibrium, e.g., by altering salt concentration and using G4-targeted ligands.

Effect of osmolytes on the thermal stability of proteins: replica exchange simulations of Trp-cage in urea and betaine solutions.

Although osmolytes are known to modulate the folding equilibrium, the molecular mechanism of their effect on thermal denaturation of proteins is still poorly understood. Here, we simulated the thermal denaturation of a small model protein (Trp-cage) in the presence of denaturing (urea) and stabilizing (betaine) osmolytes, using the all-atom replica exchange molecular dynamics simulations. We found that urea destabilizes Trp-cage by enthalpically-driven association with the protein, acting synergistically with temperature to induce unfolding. In contrast, betaine is sterically excluded from the protein surface thereby exerting entropic depletion forces that contribute to the stabilization of the native state. In fact, we find that while at low temperatures betaine slightly increases the folding free energy of Trp-cage by promoting another near-native conformation, it protects the protein against temperature-induced denaturation. This, in turn, can be attributed to enhanced exclusion of betaine at higher temperatures that arises from less attractive interactions with the protein surface.

Molecular basis of the osmolyte effect on protein stability: a lesson from the mechanical unfolding of lysozyme.

Osmolytes are a class of small organic molecules that shift the protein folding equilibrium. For this reason, they are accumulated by organisms under environmental stress and find applications in biotechnology where proteins need to be stabilized or dissolved. However, despite years of research, debate continues over the exact mechanisms underpinning the stabilizing and denaturing effect of osmolytes. Here, we simulated the mechanical denaturation of lysozyme in different solvent conditions to study the molecular mechanism by which two biologically relevant osmolytes, denaturing (urea) and stabilizing (betaine), affect the folding equilibrium. We found that urea interacts favorably with all types of residues via both hydrogen bonds and dispersion forces, and therefore accumulates in a diffuse solvation shell around the protein. This not only provides an enthalpic stabilization of the unfolded state, but also weakens the hydrophobic effect, as hydrophobic forces promote the association of urea with nonpolar residues, facilitating the unfolding. In contrast, we observed that betaine is excluded from the protein backbone and nonpolar side chains, but is accumulated near the basic residues, yielding a nonuniform distribution of betaine molecules at the protein surface. Spatially resolved solvent-protein interaction energies further suggested that betaine behaves in a ligand- rather than solvent-like manner and its exclusion from the protein surface arises mostly from the scarcity of favorable binding sites. Finally, we found that, in the presence of betaine, the reduced ability of water molecules to solvate the protein results in an additional enthalpic contribution to the betaine-induced stabilization.

Molecular dynamics simulations reveal the balance of forces governing the formation of a guanine tetrad-a common structural unit of G-quadruplex DNA.

G-quadruplexes (G4) are nucleic acid conformations of guanine-rich sequences, in which guanines are arranged in the square-planar G-tetrads, stacked on one another. G4 motifs formin vivoand are implicated in regulation of such processes as gene expression and chromosome maintenance. The structure and stability of various G4 topologies were determined experimentally; however, the driving forces for their formation are not fully understood at the molecular level. Here, we used all-atom molecular dynamics to probe the microscopic origin of the G4 motif stability. By computing the free energy profiles governing the dissociation of the 3'-terminal G-tetrad in the telomeric parallel-stranded G4, we examined the thermodynamic and kinetic stability of a single G-tetrad, as a common structural unit of G4 DNA. Our results indicate that the energetics of guanine association alone does not explain the overall stability of the G-tetrad and that interactions involving sugar-phosphate backbone, in particular, the constrained minimization of the phosphate-phosphate repulsion energy, are crucial in providing the observed enthalpic stabilization. This enthalpic gain is largely compensated by the unfavorable entropy change due to guanine association and optimization of the backbone topology.